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Recent high-resolution simulations together with theoretical studies of the dy- 
namical evolution of galactic disks have shown that contrary to wide-held be- 
liefs a 'live', dynamically responsive, dark halo surrounding a disk does not 
stabilize the disk against dynamical instabilities. We generalize Toomre's Q 
stability parameter for a disk-halo system and show that if a disk, which would 
be otherwise stable, is embedded in a halo, which is too massive and cold, the 
combined disk-halo system can become locally Jeans unstable. The good news 
is, on the other hand, that this will not happen in real dark haloes, which are 
in radial hydrostatic equilibrium. Even very low-mass disks are not prone to 
such dynamical instabilities. 
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1. Introduction 

The classical paradigm is that at given total mass of the system a galactic 
disk is stabilized against local dynamical Jeans instabilities, if it is embed- 
ded in a dark halo. This can be seen for instance from Toomre's Qt stability 
index for a stellar disk, 

3.36 GSd' ^' 
where k denotes the epicyclic frequency of the orbits of the stars, Cd is the 
radial velocity dispersion of the stars, and Ed the surface density of the 
disk. G denotes the constant of gravitation. If all other parameters are kept 
constant, but Sd is lowered, Qt rises and the disk becomes more stable 
against Jeans instabilities. The physical reasoning is that the self-gravity 
of the disk, which has a destabilizing effect, is reduced by the surrounding 
halo. Similarly the onset of non-axisymmetric coherent large-scale instabil- 
ities of the entire disk such as the bar instability was thought to be damped 
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by a surrounding halo. Ostriker & Pccblcs^^ showed in their classical nu- 
merical simulations of the dynamical evolution of a self-gravitating disk 
that the bar instability could be suppressed, if the disk was embedded in a 
halo potential. However, modern high-resolution simulations in which the 
surrounding halo is treated as a dynamically responsive system, have shown 
that actually the opposite is true. Athanassoula^ showed that in her simu- 
lations of the bar instability the bar grows stronger, if the disk is embedded 
in a live dark halo rather than in a static halo potential. This was explained 
there and, particularly, in Athanassoula^ as due to the effect of the live halo 
on the angular momentum exchange within the galaxy. First doubts about 
an entirely passive role of the halo were already raised by Toomre^^ . These 
findings were supported by theoretical studies of the swing amplification 
of shearing spiral density waves which is also enhanced, if the disk is em- 
bedded in a live dark halo instead of a static potential (Fuchs^ , Fuchs & 
Athanassoida^") . 

Here we return to the severe local Jeans instability of a self-gravitating 
disk and investigate the effect of the presence of a live dark halo. In the 
next section we demonstrate how Toomre's concept of the Qt parameter 
can be generalized in order to take into account the effect of such a halo. 
In the final section we discuss implications for realistic disk-halo systems 
(cf. also Esquivel & Fuchs''). 

2. Modification of the Qt stability index 

We study the Jeans instability of an infinitesimally thin galactic disk using 
the model of a patch of the galactic disk developed by Toomre^'' , Gol- 
dreich & Lynden-Bell-'^^ , and Julian & Toomre^^ (cf. also Fuchs^). The 
patch is assumed to rotate around the galactic center and the differential 
rotation of the stars is approximated as a linear shear flow. The surface 
density is assumed to be constant over the patch. Polar coordinates are ap- 
proximated by pseudo Cartesian coordinates {x, y) with x pointing in the 
radial direction and y in the direction of rotation, respectively. Toomre^^ 
(cf. also Lin^^) has calculated the dynamical response of the disk to a small 
'ring-like' perturbation of the gravitational disk potential of the form 

by solving the linearized Boltzmann equation which describes the evolu- 
tion of the phase space distribution of the disk stars. The induced density 
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perturbation can be expressed as 
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i{(jjt-\-kx) 



In 



OO 



1 — exp 



(3) 



2(1)^ 



$ke 



where the !„ denote modified Bessel functions and Sd the background sur- 
face density of the disk, respectively. In deriving eq. (3) a Schwarzschild 

velocity distribution of the stars with the radial velocity dispersion has 
been adopted. The disk is assumed to be self-gravitating, so that the density 
- potential pair has to fulfill the Poisson equation implying 

2wG, 



\k\ 



(4) 



Equations (3) and (4) can be combined to a dispersion relation u){k) for the 
'ring'-like perturbations of the disk. In the limit w = eq. (3) reduces to 



1 — exp 



•lo 



ikx 



(5) 



This defines a line in a space spanned by Qt and the wavelength of the 
perturbation A = 27r/fc expressed in units of Acrit = ^n^GY^d/ k'^ which, as 
can be seen by a Taylor expansion of eq. (3), separates neutrally stable 
(w^ > 0) from exponentially unstable (cj^ < 0) perturbations. The crite- 
rion that ensures that all perturbations are neutrally stable is the famous 
Toomre criterion 



Qt > 1. 



(6) 



The model of a local patch of a galactic disk has been extended by Fuchs^ 
by embedding it into a dark halo. All density gradients in the halo are 
neglected as in the disk so that the halo density distribution is assumed to 
be homogeneous. The dark matter particles follow straight-line orbits with 
an isotropic velocity distribution modelled by a Gaussian distribution. We 
can directly apply the results of Fuchs^ . As is shown there the dark matter 
halo responds to the potential perturbation 4>d in the disk and develops 
potential perturbations $h which have the same radial structure exp(jA:a;) 
as in the disk. Equations (26) and (28) of Fuchs" have been obtained as 
solutions of the linearized Boltzmann equation for the dark matter particles 
and give the Fourier coefficients of the potential perturbation in the halo at 
the midplane of the disk. In the case /cy = 0, A:|| = k, and w = 0, which we 
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consider here, there is only a contribution from non-resonant dark matter 
particles, 



where ph and cth denote the density of the dark halo and the velocity dis- 
persion of the dark matter particles, respectively. 

This induced perturbation of the gravitational potential of the dark halo 
has to be taken into account on the rhs of eq. (5), 



which means that the halo supports the perturbation of the disk and the 
density perturbation in the disk is stronger than in an isolated disk. Com- 
bining eqns. (5), (7), and (8) leads to an implicit equation that describes 
the line of neutrally stable {lo = 0) perturbations in the space spanned 
by Qt and X/Xcnt as in the case of an isolated disk, but now modified by 
the extra halo term according to eqns. (7) and (8). This can be cast into 
dimensionless form as 



with the parameters a = (3.36/27r)2 = 0.286 and /3 = {2'KGph/ul) x 
(27rG'Erf)V«-A 



In Fig. 1 we illustrate solutions of eq. (9) for various values of /3. The case 
/3 = reproduces Toomre's classical result. The unstable perturbations 
(w^ < 0) are located in the parameter space below the line. Thus for Qt > 1 
all perturbations are neutrally stable (w^ > 0). This is no longer the case, 
if finite values of /? are considered. The graphs of the solutions shown in 
Fig. 1 always turn upwards at large wavelengths A. Thus at large enough 
wavelengths'^ all perturbations of the disk - halo system become unstable 
and grow exponentially. This behaviour is related to the Jeans collapse of 
the halo component. Its Jeans length is given by Aj = ^/ira'^/Gph. or 





^ *dk + ^hk oc $dk 



(8) 





1 



(10) 



''The roots of eq. (9) are A = and, if /3 < 0.25, A/Acrit = (1±V1 - 4/3)/2/3, respectively. 
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Fig. 1. Sepciration of stable from unstable perturbations of a self-gravitating disk em- 
bedded in a live dark halo. Qt denotes the usual Toomrc stability index and A is the 
wavelength of the perturbation measured in units of A^rit. Unstable perturbations are 

located in the parameter space below the dividing lines. Lines arc shown for values of 
the /3-parameter, which describes the dynamical responsiveness of the dark halo, /3 = 0, 
0.2, 0.25, and 0.3, respectively. 



In real haloes the Jeans length will be of the order of the size of the halo or 
even larger, because otherwise the haloes would have collapsed to smaller 
sizes. Thus the inferred instability of the disk - halo system on large scales 
seems not to occur in real galaxies. As can be seen from Fig. 1 stability on 
small scales can be ensured by Qt indices at thresholds which are slightly 
larger than in isolated disks. 
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3. Discussion and Conclusions 

As a first application of the stability criterion derived here we test the sta- 
bility of the Milky Way disk and the surrounding dark halo in the vicinity 
of the Sun. The local disk and halo parameters listed in Table 1 imply 
Qt = 2.8 and (3 = 0.0057, respectively. If we include in our estimate the 
cold interstellar gas with a local surface density of 4Mq/pc^ (Dame'"') and 
assume a velocity dispersion of the interstellar gas of ag = 5 km/s, which 
leads to a reduced mass weighted effective velocity dispersion of the com- 
bined stellar and gaseous disks, the parameter values change to Qt = 2.2 
and /3 = 0.0078, respectively. Equation (10) implies that Aj = SAcrit = 39 
kpc. Since the instability sets in nominally at A = 127Acrit the Milky Way 
disk and halo system seems to be very stable. 



Table 1. Local parameters of the Milky Way 





38 Mq/pc^ 


(Holmberg & Flynn-'^^) 


o-d 


40 km/s 


(Jalircifi & Wielen") 


K 


V2- 220 km/s/8.5 kpc 


(flat rotation curve) 


Ph 


0.01 Mq/pc^ 


(Bahcall & Soneira^) 


o-h 


220 km/s/V2 


(isothermal sphere) 


-^crit 


4.8 kpc 




Aj 


39 kpc 


(/3 = 0.0078) 



In order to explore in what range the /3-parametcr of spiral galaxies 
is to be expected in general, we consider the model of a Mestel disk with 
the surface density = So embedded in a singular isothermal sphere 
representing the dark halo with the density distribution ph = poR~'^. The 
rotation curve of the model galaxy is given by 

v',{R) = vl{R) + vl{R) (11) 

with the disk contribution v'^{R) = 2TrGT,dR = const, and the halo con- 
tribution v'^{R) = AirGpiiR^ = const. (Binney & Tremaine^). Assuming 
isotropic velocity dispersions of the dark matter particles it follows from 
the radial Jeans equation, 

1 dpi,al ^ _vl 
Ph dR R ' 

that their velocity dispersion is given by 

2 -'- / 2 I 2\ 
= +^^h) , 



(12) 

(13) 
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because the dark matter partielcs are bound by both the gravitational disk 
and halo potentials. We find then 

^ R^vl + vl A {vl + vlf 4 {vl + vlf ' ^ ' 

which implies the maximal value 

/3</3ma.(i^d = 2i^h)= 0.037. (15) 

This means that in realistic halo models its density cannot be increased, 
on one hand, and the velocity dispersion of the halo particles lowered, on 
the other hand, indiscriminately because the halo has to stay in radial 
hydrostatic equilibrium. Equation (15) implies Aj = 3.7Acrit, whereas the 
instability sets nominally in at A = 26 Acrit • In order to ensure stability 
at smaller wave lengths the Toomre stability index must be larger than 
Qt > 1.02. 

As a further concrete example we analyse the dynamics of the low sur- 
face brightness galaxy F568-1. De Blok et al.^ have observed its rotation 
curve and constructed several concurring mass models for the galaxy, which 
all fit the rotation curve equally well. Using the parameters of the models 
with isothermal haloes and with either a 'realistic' M/L-ratio of the disk 
or a 'maximum- disk' we have solved numerically the radial Jeans equation 
(12) for the dark matter particles, 

2 



'lim'-\s-J'"^R'-^^ (16) 

Ph(-K) Jr R 



The resulting /3-parameters of both models are shown in Fig. 2. As can be 
seen from Fig. 2 the /3-parameters are of the same small order of magnitude 
as in the Milky Way or the simple disk - halo model described above. 

We conclude from this discussion that embedded galactic disks are not 
prone to Jeans instabilities, provided their Toomre stability index is a few 
percent higher than Qt = 1- From a practical point of view the destabilizing 
effect of the surrounding dark halo on the Jeans instability of the embedded 
galactic disks seems to be negligible. 
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Fig. 2. yS-parameters of mass models of the low surface brightness galaxy F568-1. 
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